The Flattened, Rotating Molecular Gas Core of Protostellar Jet HH 212 
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T— I ■ 

\ The recently discovered protostellar jet known as HH212 is beautifully symmetric, 

I with a series of paired shock knots and bow shocks on either side of the exciting source 

region, IRAS 05413-0104 (Zinnecker et al. 1998). We present VLA ammonia maps of 
■ the IRAS 05413-0104 molecular gas envelope in which the protostellar jet source is em- 

^ ■ bedded. We find that the envelope, with mass of 0.2 Mq detected by the interferometer, 

. is flattened perpendicular to the jet axis with a FWHM diameter of 12000 AU and an 

Q^, axis ratio of 2:1, as seen in NII3 (1,1) emission. There is a velocity gradient of about 

^ I 4-5 km sec~^ pc~^ across the flattened disk-like core, suggestive of rotation around an 

|H I axis aligned with the jet. Flux- weighted mean velocities increase smoothly with radius 

^ ' with a roughly constant velocity gradient. In young (Class 0) systems such as HH212, 

^ • a significant amount of material is still distributed in a large surrounding envelope, and 

thus the observable kinematics of the system may reflect the less centrally condensed, 
youthful state of the source and obscuration of central dynamics. The angular momen- 
tum of this envelope material may be released from infalling gas through rotation in the 
HH212 jet, as recent observations suggest (Davis et al. 2000). A blue-shifted wisp or 
bowl of emitting gas appears to be swept up along the blue side of the outflow, possibly 
lining the cavity of a wider angle wind around the more collimated shock jet axis. Our 
ammonia (2,2)/(l,l) ratio map indicates that this very cold core is heated to 14 Kelvin 
degrees in a centrally condensed area surrounding the jet source. This edge-on core and 
jet system appears to be young and deeply embedded. This environment, however, is 
apparently not disrupting the pristine symmetry and collimation of the jet. 

Subject headings: star: formation — stars: pre-main sequence — ISM: Herbig-Haro 
objects — ISM: jets and outflows — radio lines: ISM 
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1. Introduction 

The youngest protostars are often so deeply embedded in their parental cloud cores that 
shocks driven by their jets cannot be easily traced in visible shock line emission but are bright in 
near-infrared molecular hydrogen shock lines. Such is the case for the recently discovered highly 
collimated jet HH 212, near the Horsehead Nebula region of Orion (Zinnecker, McCaughrean, &: 
Rayner 1998). 

HH 212 is remarkable because of its high coUimation and symmetric pairs of bow shocks and 
shock knots on either side of the driving source (Figure 1). Theories of protostcllar jet structure 
often involve a heavy influence of the surrounding medium, with impacts creating shocked globules 
along the jet, sometimes even redirecting it. But the symmetry of HH 212 suggests that some 
pulsing mechanism from the source itself is responsible for the nearly periodic bow shock pairs. 
Somehow the symmetry is maintained even though the jet is moving through a dense gas core. In 
the innermost region, water masers have been found in motion along the outflow (Claussen et al. 
1998). 

The inferred position of the jet exciting source coincides with IRAS point source 05413-0104, 
which was detected at 25, 60, and 100 microns; it is still too cold and embedded to be detected 
strongly at 12 microns (see Beichman et al. 1986). It is surrounded by a cold (11.5 K), compact 
ammonia core out of which it presumably formed (Wouterloot, Walmsley, &; Henkel 1988). The 
spectral energy distribution climbs to longer wavelengths. It is detected also at 1.3 mm, with 
a submillimeter to bolometric luminosity of about 2% (Chini et al. 1997). It has narrow NH3 
linewidths compared to emission found at the position of other maser sources. These characteristics 
suggest that the source is very young. 

HH212 is classified as a Class system, the youngest class of protostars with much material to 
be accreted still in the surrounding molecular gas envelope (Andre, Ward-Thompson, & Barsony 
2000). Class sources have a high rate of submillimeter to bolometric luminosity (over 0.5%) and 
high rates of infall and outflow (Andre, Ward-Thompson, &: Barsony 1993, Adams, Lada, &: Shu 
1987). The morphology and kinematics of the surrounding gas at this early stage of the formation of 
a star are crucial to understand as an important link between prestellar gas cores and the resulting 
young T Tauri star /disk systems. Our VLA ammonia observations were conducted as a study of 
the gas density, temperature, and velocity distributions in this pristine, highly embedded young 
source. 
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2. Observations 

Observations were carried out with the Very Large Array (VLA) of the National Radio As- 
tronomy Observatory^. The ammonia (J,K) = (1,1) and (2,2) inversion transitions were observed 
simultaneously, at rest frequencies of 23.694495 GHz and 23.722733 GHz, respectively. The width 

of the antenna primary beam response was approximately 2' FWHP. The array was in its compact, 
"D" configuration. The synthesized beamsize was 8.4" x 8.8" FWHM, with natural weighting and 
a 20 kA taper applied to the UV data. The bandwidth was 3.125 MHz and the channel separation 
was 24.4 kHz, corresponding to 0.31 km sec^^. 

The absolute flux calibrator was 3C286, with a calculated flux density of 2.41 Jy at 1.3 cm. 
The bandpass calibrator was 3C84, and a phase closure solution was applied based on interleaved 
observations of the phase calibrator B0605-085 (1950). The "dirty" beam was deconvolved from 
the maps using the "Imager" algorithm in the NRAO AIPS image processing software. 

3. A Flattened, Rotating, and Heated Gas Envelope 

The integrated intensity of the NH3 (1,1) line is displayed in the contours of Figure 1, overlaid 
on the 2.12 /im. image of HH 212 from Zinnecker, McCaughrean, & Rayner (1998). The gas core is 
flattened with an axis ratio of about 2:1. The major axis has a FWHM extent of 29" , corresponding 
to 12000 AU at an assumed distance of 400 pc, with full extent of over 15000 AU; it is centered on 
the jet source, and it lies perpendicular to the jet which has a position angle of about 25 degrees 
east of north and only 4dcg from the plane of the sky (Claussen et al. 1998). 

A velocity gradient is apparent across the flattened gas core. In Figure 2a the integrated inten- 
sity (moment 0) is again plotted in contours, while the underlying shades represent the (moment 
1) velocity field. Radial velocities change smoothly from blue-shifted gas in the northwest side to 
redshifted gas in the southeast; the large velocity gradient of about 4-5 km sec~^pc~^is virtually 
constant. The direction of the gradient within the central FWHM intensity region of the core is 
perpendicular to the jet axis, with the nodes of constant velocity roughly parallel to the jet. 

Along the eastern side of the core, a small amount of blueshifted material is seen to extend to 
the northeast alongside the northern, blueshifted side of the jet. To a lesser extent, the same effect 
is also seen to the south with a small protrusion of redshifted gas extending in the same direction 
as the redshifted side of the jet. Molecular gas is possibly being swept along by the outflow. There 
is even a hint, at least in the north, that molecular gas is found in a "bowl" surrounding the jet 
axis, possibly outlining a cavity of a wider-angle flow surrounding the coUimated shock emission of 
the jet (see, e.g., Shu et al. 1995, Shang, Shu,&; Glassgold 1998). 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 
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Ammonia inversion lines are a convenient probe of both optical depth and temperature. Most 

ammonia inversion transitions are composed of multiple hyperfine components. The relative inten- 
sities of central and satellite component blends of the (1,1) transition can be used to calculate the 
optical depth. For the main component and total transition of the (1,1) line, we find r values of 0.8 
and 1.6, respectively, at the core center. Using the procedures and assumptions described in Ho k, 
Townes 1983 (see also Wiseman & Ho 1998), we also estimate an NH3 column density of 1.6 x 10^^ 
cm~^ and a total volume density n of 1 x 10^ cm~^. The mass of the core gas detected with the 
interferometer is approximately 0.2 M©, assuming an NH3/H2 abundance ratio of 2 x 10~® (Harju, 
Walmsley, & Wouterloot 1993). Within the large uncertainties of this abundance ratio, the mass 
is the same as that needed to gravitationally bind the rotating system. 

The ratio of intensity of the NH3 (2,2) and (1,1) transitions reflects level populations and thus 
serves as a thermometer for rotational temperature within cold molecular gas cores. In Figure 2b, 
the contours again show the integrated intensity of (1,1) emission. The shades show the (2,2)/(l,l) 
ratio, indicating the relatively warmer regions. It is apparent that this is indeed a cold core of gas, 
with the (2,2) emission showing up only in the central region of the core. The (2,2)/(l,l) ratio 
peaks near the position of the embedded protostellar source, and extends slightly to the northwest. 
The peak rotational temperature here is 14 K, as calculated from the line ratio (Ho & Townes 
1983). In colder regimes such as this, the kinetic temperature is approximately equal to the NH3 
rotational temperature (Danby et al. 1988, Kuiper 1994). The young protostar is evidently still in 
the center of the gas core, and we may be seeing here the warm gas heated by the protostar and 
the jet as they interact with the environment. 



4. Discussion 

4.1. Rotation in Protostellar Systems 

Envelopes of Class protostars carry much (often most) of the mass of the system and can 
extend to diameters of several thousands of AU. Interferometric maps of such systems tend to 
show flattened disks perpendicular to the axis of the associated jet (Myers, Evans, k, Ohashi 2000). 
The flattened molecular gas envelope of the HH 212 jet soTuce has a smooth and roughly constant 
velocity gradient perpendicular to the jet axis (Figure 2a, 2c). 

Turbulence in cloud cores can produce velocity gradients that mimic rotation (Burkcrt & 
Bodenheimer 2000 (BB)). However, the gradient we find in HH 212 is faster than the most likely 
effect produced in the BB turbulence models. It should be noted that the BB models were based on 
more massive cores, so proper scaling of the models would be required to better assess the possibility 
of a turbulent component to the observed velocity gradient in HH 212. But it is also true that the 
HH 212 gradient extends systematically across the long axis of the core, and perpendicular to the 
jet, for over 5 synthesized beamwidths. Thus the combined evidence suggests that the gradient is 
more likely dominated by rotation. 



-5- 



The more evolved T Tauri systems often show Keplerian rotation signatures in the less plentiful 
envelope gas, with the radial velocities increasing and diverging toward the central source {V oc 
An early report of Keplerian rotation in a young stellar disk was made for -"^^CO observations 
of HL Tau (Sargent &: Beckwith 1987, 1991), though in this source the kinematics proved to be 
much more complicated, with evidence for infall (Hayashi et al. 1993) and outflow (Cabrit et al. 
1996). The strongest Keplerian velocity curves are seen in circumstellar disks of a few hundred 
AU in Class II T Tauri systems without much envelope material left, thus having little confusion 
from large scale infall or outflow motions (Mundy, Looney, & Welch 2000, Dutrey et al. 1998). 
Keplerian rotation has been detected in the gas disks of GM Aur (Dutrey et al. 1998) and DM Tau 
(Guilloteau & Dutrey 1998), for example, as well as in multiple systems such as UY Aur and the 
beautiful ring of GO Tau (Duvert et al. 1998, Guilloteau, Dutrey, &; Simon 1999). 

Several Class systems such as HH 212 are in dense gas cores that appear to have a more 
rigid rotation signature (Myers, Evans, Sz Ohashi 2000); this is consistent with their character 
of having much of the system mass still distributed in the gas envelope around the embedded 
accreting protostar(s). Indeed, the 0.2 Mq we detect with the VLA in the HH 212 core is 10 times 
the circumstellar gas mass found in some T Tauri systems (e.g. < O.OIMq for UY Aur (Duvert 
et al. 1998)). In such early systems, a gradient in the intensity-weighted mean velocity that is 
constant rather than diverging toward the center can indicate that either a) the non-Keplerian 
rotation curve is reflecting a mass distribution that is not yet dominated by a central condensation, 
or b) the system is in Keplerian rotation, but high optical depths result in velocity maps that are 
dominated by motion in the outer layers along the line of sight toward the center of the core. Models 
of expected channel map emission can be used to test for Keplerian rotation (e.g. Kocrncr, Sargent, 
& Beckwith 1993, Dutrey, Guilloteau, & Simon 1994). Higher spatial and spectral resolution maps 
would be needed to fully model the HH 212 envelope, although the slightly broader line widths we 
observe toward the center of the core are possible evidence for some Keplerian component to the 
motion. In any case, it now appears that the observed kinematics of the protostellar envelope can 
signify the developmental stage of the embedded YSO. More high resolution kinematical studies are 
needed to confirm this pattern. The large envelopes such as we observe in HH 212 may feed mass 
directly to the embedded protostar (if at the earliest stage) or to the inner circumstellar accretion 
disk just beginning to form. It is important to note that this inner disk (< 100 AU and thus 
unobservable in our NH3 maps) probably rotates with near-Keplerian velocities (Stahler 2000). 

Many rotating cores of these large sizes (> 12000 AU) do not appear to conserve specific 
angular momentum, defined as Vrot x Rrot) where Vrot is the velocity of rotation at radius Rrot 
(Myers, Evans, & Ohashi 2000). Our results show that, with a specific angular momentum of 
0.012 km s~^ pc at a radius of 6000 AU, the HH 212 gas envelope would fit well with other cores 
of similar size in a plot of Ohashi et al. (1997) showing specific angular momentum falling with 
decreasing source radius. However, Ohashi ct al. showed that there is a transition at about this 
same size scale; smaller dense core regions do show conserved specific angular momentum in regions 
of both massive infall and rotationally supported disks. Thus this would serve as the size scale for 
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dynamical collapse, where angular momentum is conserved. The HH 212 core seems to be in fast, 
bound rotation at the size scale of this transition. 

Since the velocity gradient and the long axis of the envelope are both perpendicular to the 
jet axis, it is suggestive that rotation may contribute to the observed flattening. We calculate 
rotational energy of the system to reach up to about 3-10% of each of the thermal, turbulent, 
and gravitational energy (Harju et al. 1993, Myers 1983). Preferential contraction along an axis 
perpendicular to the disk plane, along possible magnetic field lines, could also produce flattened 
morphology. Another intriguing possibility is that the outflow itself contributes to the morphology 
by clearing out dense gas along the flow. Indeed, the blue-shifted wisp of gas extending to the 
northeast alongside the (northward) blue-shifted side of the jet (Figure 2a) seems to form part of 
a "bowl" around the jet and the associated CO outflow recently studied by Lee et al. (2000); part 
of the blue lobe of the CO flow shifts to the east corresponding to the NH3 wisp. It has been 
suggested that molecular outflows associated with jets may broaden over time, as has been seen 
in B5-IRS1; the broadening flow could contribute to the eventual decline of accretion (Velusamy 
k. Langer 1998). Thus we may be seeing the HH 212 flow in an early stage of broadening and 
sweeping out envelope material. 



4.2. Angular Momentum Release Through the Jet 

Bipolar outflows appear to be a ubiquitous product of the formation of stars. Strong winds 
from young stellar sources are collimated into jets which are observed in visible and infrared shock 
emission. Magnetic fields may be responsible for channeling and collimating the jets (see, e.g., Shu 
et al. 1995). 

Outflows are thought to be a 'release valve' for the angular momentum carried into the proto- 
stellar accretion region by infalling material. Otherwise, matter infalling from even a very slowly 
rotating circumstellar envelope will have too much conserved angular momentum to fall into the 
inner accretion region, unless sufficient angular momentum is lost through disk viscosity. 

Observational evidence for jets as significant angular momentum carriers has been virtually 
nonexistent. Recently, however, results of slit spectroscopy on the shock knots in the HH 212 jet 
have given some evidence of jet rotation (Davis et al. 2000). Five of six knots studied show similar 
evidence of rotation; in the inner part of the southern flow, for example, a clear velocity gradient 
is seen. If this gradient represents rotation, the jet rotation is in the same directional sense as 
that of the rotating molecular gas envelope we observe in ammonia. A rough extrapolation of the 
large scale angular momentum in the envelope gas fits the observed gradient in the jet, assuming 
a fraction of infalling mass is channeled into the fiow. Thus this may be the first observational 
evidence for a rotating jet carrying away the angular momentum of material falling in from a 
rotating envelope. Clearly more detailed observations and evidence of jet rotation and inner disk 
kinematics are needed. 



- 7- 



We thank our referee and also P. T. P. Ho for insightful guidance. Support for J.W. was 
provided by NASA through Hubble Fellowship grant 7^HF-01115.01-98A awarded by the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in 
Astronomy, Inc., for NASA under contract NAS 5-26555. 

REFERENCES 

Adams, F. C, Lada, C. J. & Shu, F. H. 1987, ApJ, 312, 788 
Andre, P. , Ward-Thompson, D. k Barsony, M. 1993, ApJ, 406, 122 

Andre, P., Ward-Thompson, D. & Barsony, M. 2000, Protostars and Planets IV (Book - Tucson: 
University of Arizona Press; eds Mannings, V., Boss, A. P., Russell, S. S.), p. 59, 59 

Beichman, C. A., Myers, P. C, Emerson, J. P., Harris, S., Mathieu, R., Benson, P. J. & Jennings, 
R. E. 1986, ApJ, 307, 337 

Burkert, A. & Bodenheimer, P. 2000, ApJ, 543, 822 

Cabrit, S., Guilloteau, S., Andre, P., Bertout, C, Montmerle, T. and Schuster, K. 1996, A&A, 305, 
527 

Chini, R., Reipurth, B., Sievers, A., Ward-Thompson, D., Haslam, C. G. T., Kreysa, E. and Lemke, 
R. 1997, A&A, 325, 542 

Claussen, M. J., Marvel, K. B., Wootten, A. & Wilking, B. A. 1998, ApJ, 507, L79 

Danby, G., Flower, D. R., Valiron, P., Schilke, P. and Walmsley, C. M. 1988, MNRAS, 235, 229 

Davis, C. J., Berndsen, A., Smith, M. D., Chrysostomou, A. and Hobson, J. 2000, MNRAS, 314, 
241 

Dutrey, A., Guilloteau, S., Prato, L., Simon, M., Duvert, G., Schuster, K. &; Menard, F. 1998, 
A&A, 338, L63 

Dutrey, A., Guilloteau, S., & Simon, M. 1994, A&A, 286, 149 

Duvert, G., Dutrey, A., Guilloteau, S., Menard, F., Schuster, K., Prato, L. h Simon, M. 1998, 
A&A, 332, 867 

Guilloteau, S. & Dutrey, A. 1998, A&A, 339, 467 

Guilloteau, S., Dutrey, A. & Simon, M. 1999, A&A, 348, 570 

Harju, J., Walmsley, C. M. & Wouterloot, J. G. A. 1993, A&AS, 98, 51 

Hayashi, M., Ohashi, N. & Miyama, S. M. 1993, ApJ, 418, L71 



-8- 



Ho, P. T. P. and Townes, C. H. 1983, ARA&A, 21, 239 

Koerner, D. W., Sargent, A. I., &: Beckwith, S. V. W. 1993, Icarus, 106, 2 

Kuiper, T. B. H. 1994, ApJ, 433, 712 

Lee, C, Mundy, L. G., Reipurth, B., Ostriker, E. C. & Stone, J. M. 2000, ApJ, 542, 925 

Mundy, L. G., Looncy, L. W., & Welch, W. J. 2000, in Protostars and Planets IV, eds. V. Mannings, 
A. P. Boss, & S. S. Russell (Tucson: University of Arizona Press), 355. 

Myers, P. C. 1983, ApJ, 270, 105 

Myers, P. C, Evans, N. J., Sz Ohashi, N. 2000, in Protostars and Planets IV, eds. V. Mannings, 
A. P. Boss, & S. S. Russell (Tucson: University of Arizona Press), 217. 

Ohashi, N., Hayashi, M., Ho, P. T. P., Momose, M., Tamura, M., Hirano, N. and Sargent, A. I. 
1997, ApJ, 488, 317 

Sargent, A. I. & Beckwith, S. 1987, ApJ, 323, 294 

Sargent, A. I. k Beckwith, S. V. W. 1991, ApJ, 382, L31 

Shang, H., Shu, F. H. and Glassgold, A. E. 1998, ApJ, 493, L91 

Shu, F. H., Najita, J. , Ostriker, E. C. & Shang, H. 1995, ApJ, 455, L155 

Stahler, S. W. 2000, in "Star Formation from the Small Scale to the Large Scale," eds. F. Favata, 
A. A. Kaas & A. Wilson (ESA: SP-44), p. 133. 

Velusamy, T. & Langer, W. D. 1998, Nature, 392, 685 

Wiseman, J. J. and Ho, P. T. P. 1998, ApJ, 502, 676 

Wouterloot, J. G. A., Walmsley, C. M. and Henkel, C. 1988, A&A, 203, 367 
Zinnecker, H., McCaughrean, M. J. & Rayner, J. T. 1998, Nature, 394, 862 



This preprint was prepared with the AAS IM^^X macros v5.0. 



-9- 



Fig. 1. — The HH 212 jet, as imaged in H2 2.12 micron shock emission (Zinnecker, McCaughrean, 
& Rayner 1998). The observations were carried out in Spain in 1997 with the Calar Alto 3.5-m 
telescope using its Omega Prime wide field infrared camera and kindly given courtesy of T. Stanke. 
The contours represent the integrated NH3 (1,1) emission of the dense protostellar envelope gas, 
as mapped with the VLA. Contour levels begin at 16 mJy beam~^ km s~^ and increase by steps 
of 16 mJy beam~^ km s~^. 

Fig. 2. — (a): The velocity field of the dense gas around the HH 212 jet, as imaged with the VLA 
in the ammonia (1,1) line. The shades represent the moment 1 velocity field, showing a velocity 
gradient across the flattened gas envelope as evidence of rotation. Contours are the integrated 
intensity as in Figure 1. The arrows indicate the orientation of the HH 212 jet. (b): Heating 
in the HH 212 NH3 core. Contours arc the integrated intensity of the (1,1) line, as in Figure 1. 
The shades represent the (2,2)/(l,l) intensity ratio, which increases with temperature and optical 
depth. The inner region is slightly heated to 15 K, possibly where the embedded YSO and jet are 
heating the surrounding gas. (c): Velocity-position cut across the FWHM long axis of the HH 212 
flattened envelope, showing the gradient in first moment velocity. Positions are taken across the 
midline of the flattened envelope of panel (a), with offsets taken from the core center along an axis 
perpendicular to the jet. 
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